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1. INTRODUCTION 

Spectral analysis of ultrasound Doppler signals is an important tool in the diagnosis of vascular 
diseases and in the study of stenosis. The most widely used method to analyze a Doppler signal in real time is 
the STFT short-time Fourier transform. Unfortunately, this method has a major disadvantage, which is the 
absence of a compromise between the resolution of the frequency and the resolution time. However, a good 
time and frequency resolution is necessary to characterize the non-stationarity of the signal. In addition, 
it must be assumed that the signal can be processed as a stationary stochastic process in narrower 
analysis windows. 

The basic idea of time-frequency analysis is to understand and describe situations where the frequency 
content of a signal changes over time. Although the time-frequency analysis is almost 70 years old, significant 
progress has been made over the last 20 years. In particular, time-frequency representation has received 
considerable attention as a powerful tool of high resolution and accuracy for the analysis of a variety of 
biomedical signals such as ultrasound Doppler signals, as well as signals from other fields [1], [2]. 

In this paper, we have limited our study to a few distributions such as the Wigner- Ville distribution 
(WVD), the Pseudo Wigner-Ville-distribution (PWVD), and the smoothed pseudo Wigner-Ville distribution 
(SPWVD). Other techniques based on the wavelet transform such as the S-transform can also be applied to 
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analyze the Doppler blood flow signal. The WVD, PWVD and SPWVD distributions were introduced for non 
stationary spectral analysis and showed good time-frequency resolution [3]. 


2. METHODS OF ANALYSIS 
2.1. The Fourier Transform 

Ultrasound Doppler signals are non-stationary signals, because the frequency that represents the speed 
of red blood cells varies with time, in addition Doppler signals is frequency-modulated by the speed of red 
blood cells movement. 

To perform a spectral analysis of a Doppler signal, we use the Fourier transform given by the 
following relation: 


X(f) = [0 x(e2™at (1) 


This transform provides information on the frequency components of the signal studied; their 
amplitude and their phase. The major disadvantage of this application is that it causes a loss of information 
concerning their temporal location [3]-[5]. 


2.2. The S-Transform 

The S-transform of a signal x(t) representing a time-frequency analysis x(t) was developed by Stokwell 
in 1996. This transform is similar to the short-term Fourier transform, except that the amplitude and the width 
of the analysis window is variable depending on the frequency, as is the case in wavelet analysis; this transform 
is given by: 


Sirf) = [°° x(w(t — 1, f).e7?™at (2) 


Where the mother wavelet is written by: 


~f2¢2 
If| 


w(t f) = T= ek” Jk > 0 (3) 


Where ¢ is a time variable, fis the frequency, and k is a scale factor that controls the number of 
oscillations in the mother wavelet [6]. 


2.3. The discrete S Transform 
Let us consider X /kT], k = 0; 1; ...; N-1 a discrete signal, corresponding to x(t) with a sampling period 
T. The discrete Fourier transform is given by [7]: 


n 1 N-1 _j2mnk 
x |=] = — yRodx[kT]e7 8 (4) 
Where n= 0, J..., N- J. In the discrete case, the S transform is the projection of the vector defined by 
x[kT] in a non-orthogonal basis of vectors and the elements of the S transform are not independent. 
Each basis vector (of the Fourier transform) is divided into N localized vectors by an element-by-element 
product with the shifted Gaussians, such that the sum of these N localized vectors 1s equal to the original basis 


vector. By using the discrete version of the (2), and when placing f > = and t > IT, the S transform of a 
discrete signal x/kT/] 1s given by: 


m+n 
NT 


n? e@ N (5) 
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wp | _ yN-1 
S it, ~ = nye | 
and for the n = 0 voice, 1s equal to the constant defined by: 
= NS | Se 
S[iT, 0] = —YN=3 x|=] (6) 
Where i, m and n=0, /..., N-1. 


The (6) puts the constant average of the signal in the way corresponding to zero frequency. 
Thus, the existence of the reverse 1s ensured in the general case. The discrete S transform faces the problems 
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of sampling and finite length, causing the implicit periodicity in the time and frequency fields. 
The convolution operations are implicitly ' wrap-around ' causing edge effects. The reverse of the discrete S 
transform is given by: 


i2m™k 


x[kT] = SONG {ENS [Zit |}e we 7) 


2.4. The Analytical Signal 
The analytical signal is used to avoid the overlap of the frequency components in the time-frequency 
plane; this signal is obtained by applying the Hilbert transform: 


x, (t) = x(t) + jHix(t)} (8) 


Where Hx (t)} represents the Hilbert transform of x(t) 


H{x(t)} = Ivp JS at (9) 


Fe. ee 
vp being the abbreviation of Cauchy's main value [8]. 


2.5. The Wigner-Ville Transform 

This approach was developed by Wigner for use in physics, but was later applied by Ville for signal 
processing. This Wigner-Ville transform provides time-frequency decomposition without any restriction on 
temporal and frequency resolutions. It appears to be perfectly suited to the analysis of non-stationary signals 
since it does not require assumptions about the signal itself. TWV is defined by: 


WV, (tf) = 7x (t+=)x" (t-te Prd (10) 


This formula represents the energy of a signal x at time ¢ and at frequency f/ Unfortunately, the non- 
linearity of this transform has disastrous consequences manifested by the appearance of interference and 
negative energies in the time-frequency distribution of the signal energy. These phenomena are remarkable 
when there are two distinct events in time and frequency [8]. This phenomenon is unrealistic from a physical 
point of view, and suggests difficulties in the interpretation of the TWV. Therefore, and in order to extract 
reliable results, it is necessary to avoid these negative energies by post-treatment that will affect the resolution 
properties of the transform. The DWV is related to the ambiguity function A,(n, z), which is a time-frequency 
correlation measurement; by a two-dimensional Fourier transform relation: 


A,(n,t) = [7 [7x (ut4)x*(u—S) edu = fw, erent at. df (11) 


This makes it possible to consider the classical problems involving the ambiguity function (in 
particular optimal detection-estimation, adapted filtering and signal synthesis) in a dual way by working 
directly on the time-frequency plane with the VWD. 


2.6. The Pseudo Wigner-Ville Transform 

The pseudo Wigner-Ville Distribution (PWVD) is a_ short-term Wigner-Ville transform. 
An observation time window h is introduced, and this distribution smooths out in the frequency directions only 
where the interference will not be attenuated in time. 


PWV, (t, f) = (ae h(t) x (t + =) x (t — =) e 2nd (12) 


If we introduce a temporal smoothing in the equation using a frequency window g, we will have a 
smoothed pseudo Wigner- Ville distribution [9], [10]. 


2.7. The Smoothed Pseudo-Wigner- Ville Transform 

The pseudo Wigner-Ville distribution has the effect of reducing the small parasitic ripples in the 
frequency direction. However, the parasitic ripples appearing in the PWVD can be further reduced by 
smoothing the PWVD, this time in the direction of the time of the plan (¢, f), where we have a new transform, 
i.e. that of the smoothed pseudo-Wigner- Ville distribution given by the following relation [11], [12]: 
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SPWV,(t,f) = J7 h(t) J e(s —t) x(s +4)x* (s—).ds.e Pdr (13) 


3. RESULTS AND DISCUSSION 

The major theoretical advantage of the methods studied in the previous part lies in their ability to 
define more precisely the non-stationary signals in the time-frequency space without any hypothesis on the 
characteristics of the signal. Spectral analysis of non-stationary signals requires the implementation of specific 
tools to describe a temporal evolution of frequency characteristics. In addition, they can provide precise 
information on the speed of blood flow in the femoral arteries [12]. 

The Doppler blood flow signal is a non-stationary Gaussian random signal with a relatively slow 
average variable frequency. The work of Keeton and S. Kara shows that the Ultrasonic Doppler signal of the 
femoral arteries can be considered as a stationary signal on segments from /2 ms to 20 ms; the choices of the 
size of the analysis windows 1s based on this proposition [13]. 

The ultrasound Doppler signals studied in this work are real signals from St Marie's Hospital in 
Leicester, and have a sampling frequency of 10240 HZ, and duration of 4 s. The following Figure 1 shows the 
time representation of one of these signals. 
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Figure 1. (a) Temporal representation, (b) Spectrum of doppler ultrasound signal 


A time-frequency distribution is used to analyze and process non-stationary signals in the time- 
frequency domain, as well as to improve the performance of time-frequency analysis with respect to the spectral 
analysis of biomedical signals. 

Figures (2a) and (2b) represent the sonograms obtained by PWD and SPWD of the analytical Doppler 
signal, using Hamming windows of different widths. Note, however, that small parasitic ripples are much 
smaller in PWVDs than in WVDs. It can be seen that the sonograms obtained by the PWD and SPWD have a 
better frequency resolution than the results obtained by WVD. The degraded frequency resolution is the 
disadvantage of the PWVD compared to the WVD, but the PWVD always keeps the same temporal resolution, 
since the convolution operation only acts on the frequency direction. This is interesting compared to the 
spectrogram, which provides a representation in which the frequency and time resolutions are always related. 
A judicious choice of types and widths of smoothing windows is a necessary condition for having a good 
resolution on the time-frequency plan in the application of PWVD or SPWVD on ultrasonic Doppler signals 
[14]-[16]. 

The WVD retains the temporal and frequency supports of a signal, and thus allows a precise location 
both in time and in frequency. This is another fundamental property that enables to conveniently separate the 
different events on the time-frequency plane. 
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Figure 2. Femoral arterial doppler sonograms using: (a) WVD, (b) PWVD, (c) SPWVD, (d) S-transform 


Consequently, due its bilinear structure, the WVD naturally creates interference terms between the 
different components of a signal. These terms, which are in some ways useful because they carry information 
on the relationships between the different components, are generally disruptive vis-a-vis the readability of a 
time-frequency diagram. It can indeed be shown that they are characterized by an oscillatory structure, 
in contrast with the term "signal" with softer behavior [17], [18]. 

The sonograms obtained by the SPWVD show in Figure 2(c) give better performances compared to 
other time-frequency techniques. This technique has a better resolution and a minimum of interference and 
spectral spread compared to the other techniques due to the Wigner- Ville distribution. 

Figure 2(d) illustrates the sonogram obtained by the S-transform; this result shows a quantitative 
improvement in the time-frequency analysis compared to the results obtained by Wigner- Ville distributions, 
and this is due to the use of wavelets in this technique [19], [20]. 

Note, however, that small parasitic ripples are much smaller in TPWVs than in TWVs. This aspect of 
compromise between the frequency resolution and the degree of attenuation of the parasites (interference terms) 
is in fact a typical and classic problem that is encountered in most practical situations 
(multi-component signals). This is also the reason why the use of the 7TPWV always requires a judicious choice 
of the window (in terms of type and width) according to the characteristics of the signal processed and the 
information sought [21], [22]. 


4. SBI CALCULATION 

The frequency envelopes that are extracted from the sonograms obtained by different methods will be 
applied to the input of an order 6 averaging filter to obtain a better resolution at the systolic peaks [13]. 

The degree of stenosis (given in percentage %) 1s expressed as the ratio of the diameter reduced by 
the stenosis to the actual diameter of the artery. Maximum and mean frequency envelopes of the carotid arterial 
extracted from: (a) WVD, (b) PWVD, (c) SPWVD, (d) S-transform as shown in Figure 3. 

The flow upstream of the stenosis is laminar and the red blood cells flow with a speed called average 
speed "Vinoy'. At the level of the stenosis, the speed of the red blood cells increases because of the arterial 
narrowing, to maintain a constant flow. In this case, the speed is maximal, and is called "Ving". 
The effect of stenosis on the flow of blood in the arteries as shown in Figure 4. 


Indonesian J Elec Eng & Comp Sci, Vol. 14, No. 3, June 2019 : 1220 — 1227 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 OO 1225 




































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































ssl | || H ‘| | 

a | | os | | 

F oss fll | h A Boal {IN ) i | Ne 
ml | | | j mal " i ! | 
“dhe dal =A 

g ostillit| Sane i 7 

bo | | | | Hl | | al | | | fm | 
_ a vy | rh yh mw 

| hes ) | nag! - ) 


Figure 3. Maximum and mean frequency envelopes of the carotid arterial extracted from: 
(a) WVD, (b) PWVD, (c) SPWVD, (d) S-transform 


Immediately downstream of the stenosis, when the diameter increases sharply, a swirling motion 
appears and totally disorganizes the flow; this allows the red blood cells to take multiple speeds and in 
all directions. 
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Figure 4. The effect of stenosis on the flow of blood in the arteries 


The average value of these speeds gives the average speed of the flow Vmean. All velocities measured 
at the stenosisand near it, including the maximum velocity "Vinax", are represented by the Doppler spectrum as 
frequency components. In this case, we can say that the maximum frequency of the Doppler spectrum "finax"' 
measures the maximum speed of the flow "Vinax'", and the average frequency "fmean'' Measures the average speed 
Vien. 


fmax 18 calculated by the average of the frequencies at the systolic peaks. The results obtained are 
presented as follows: 
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Relationship between SB/ and the degree of stenosis: 


SB] — fmaxcimean (14) 


fmax 


Increasing the degree of stenosis increases the flow rateat the level of arterial narrowing; this is 
interpreted by the following equation: 


B= A, (ean) (15) 


Vmax 


Where Vmean is the flow velocity when the diameter of the artery is equal to A, and Vinax becomes the 
flow velocity when the diameter is equal to B. 
Then express the degree of stenosis can be expressed as a function of Vmean and Vmax: 


(A 





—F) x 100% = Sue x 100% 


A max 


Degree of stenosis = (16) 


Since the velocities Vinay and Vinean respectively represent the frequencies finax and fmean of the Doppler 
spectrum, we can express the degree of stenosis as a function of fax and fmean: 
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Figure 5. Filtered envelopes fmax and fmean of the artery femoral sonograms using: 
(a) WVD, (b) PWVD, (c) SPWVD, (d) S-transform 


This last equation measures the spectral broadening between finax and fmean With respect to finax; which 
represents the spectral broadening index SBI. 

The frequencies finax and fnean are obtained respectively by calculating the average of the different finax 
and fmean ON a Succession of systolic cycles generated by the sonograms. 
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5. CONCLUSION 

The main objective of this work is to present the performances of several non-parametric time- 
frequency techniques, namely the distributions of Wigner-Ville, pseudo Wigner-ville, smoothed pseudo 
Wigner Ville and the S-transform in the analysis of Doppler ultrasound signals. The application of these 
methods on several Doppler signal recordings of normal Doppler arteries or showing certain pathologies allows 
us to obtain sonograms to have a better time-frequency analysis. Based on the results obtained by the different 
distributions of Wigner- Ville and the S-transform, it 1s clear that the sonograms have been reached by the WVD, 
but that the use of the SPWD makes it possible to smooth these interferences; the S-transform generates 
sonograms whose spectral and temporal resolutions are better compared to the Wigner- Ville distributions. Our 
study is followed by a quantified assessment of the severity of the SB/ study. The SB/ was calculated for 
different normal and pathological cases. Results have showen that there is a relationship between SB/ and the 
degree of surface stability. 
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